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[57] ABSTRACT 

Apparatus for obtaining significant height information 
of ocean waves, or peaks of rough terrain utilizing 
means for compressing the radar signal over different 
widths of the available chirp or Doppler bandwidths, 
and means for cross-correlating one of these images 
with each of the others, where the center frequencies of 
the images have a spacing Af*> Upon plotting a fixed 
(e.g., zero) component of the cross-correlation values 
R £ (Af,o) as the spacing is increased over some empiri- 
cally determined range, such as 0 to 1 MHz, the system 
is calibrated. Thereafter to measure height with the 
system, a spacing value is selected and a cross-correla- 
tion value is determined between two intensity images 
at a selected frequency spacing, such as 0.1 MHz. The 
measured height is the slope of the cross-correlation 
value determined to the spacing value used. Both elec- 
tronic and optical radar signal data compressors and 
cross-correlations are disclosed for implementation of 
the system. 

5 Claims, 7 Drawing Figures 
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SURFACE ROUGHNESS MEASURING SYSTEM 
ORIGIN OF INVENTION 

The invention described herein was made in the per- 5 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-568 
(72 Stat. 435; U.S.C. 2457). 

BACKGROUND OF THE INVENTION 10 

This invention relates to novel apparatus for measur- 
ing heights of peaks above an arbitrary level by statisti- 
cal means not limited to roughness or minimum area, 
and more particularly to a synthetic-aperture radar 15 
system for measuring the average height of terrain 
peaks or ocean waves. 

Much work has been done on the remote measure- 
ment of ocean wave heights using nonimaging radar, 
and a review of such techniques is presented by D. E. 20 
Barrick, Remote Sensing of Sea State by Radar, Remote 
Sensing of the Troposphere , edited by V. E. Derr, Ch. 12, 
U.S. Government Printing Office, Washington, D.C., 
1972. G. T. Ruck, D. E. Barrick, and T. Kaliszewski, 
Bistatic Radar Sea State Monitoring, Battelle Columbus 25 
Laboratories, Columbus, Ohio, Rep. 1388-19, 1971, 
have analyzed the possibility of using two-frequency 
radar interferometry for measuring wave slopes, and 
Weissman, et al. have demonstrated the operation of a 
radar system utilizing this technique both in laboratory 30 
wave tank measurements and flights over the Carib- 
bean, [D. E. Weissman, Two Frequency Radar Interfer- 
ometry Applied to the Measurement of Ocean Wave 
Height, IEEE Trans, on Ant. and Prop., AP-21, 
649=656, 1973; D. E. Weissman, C. T. Swift, W. L. 35 
Jones, Jr., J. W. Johnson, W. L. Grantham, J. A. How- 
ell, J. C. Fedors, and J. J. Davis, A Dual Frequency 
Radar for Ocean Roughness Sampling, URSI Fall 
Meeting, Boulder, Colorado, 1974.] The possibility of 
measuring ocean wave heights by studying the return 40 
from a short-pulse radar has also been demonstrated by 
B. D. Zamarayev and A. I. Kalmykov, On the Possibil- 
ity of Determining the Spatial Structure of an Agitated 
Ocean Wave by Means of Radar, Izvestia, USSR Acad- 
emy of Sciences — Atmospheric and Ocean Physics, 5, 45 
64-66, 1969, and by D. L. Hammone, R. A. Manella and 
E. J. Walsh, Short Pulse Radar Used to Measure Signif- 
icant Wave Height and Surface Wind Speed, URSI Fall 
Meeting, Boulder, Colorado, 1975. 

The principal concept involved in the wave height 50 
measurement, which is common to the interferometric 
techniques utilizing radar sensors, is that the resultant 
intensity for the electric field scattered from the ocean 
wave at some arbitrary point in space changes as the 
wavelength or angle of illumination for the electromag- 55 
netic radiation varies. This variation is dependent on the 
height difference between the crest and trough of the 
wave, and a measurement of the rate of this variation 
provides an indication of the ocean wave height. This 
principle has been used before in the design and opera- 60 
tion of the tellurometer for the precise determination of 
range differences, [T. L. Wadley, Electronic Principles 
of the Tellurometer, Trans. South African Inst. Electr. 
Engrs., 49, 143-161, 1958.] The various radar techniques 
for measuring ocean wave heights differ in the handling 65 
of the radar data, which is governed by the constraints 
set by the particular radar instrument being used, the 
accuracy and spatial resolution of the measurement 
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being related to the bandwidth and the resolution of the 
radar system involved. 

An extensive theoretical and experimental treatment 
of the utilization of this concept to measure the height 
of surface features within a resolution cell in images 
obtained by a tunable dye laser has been performed by 
Jain, the present inventor [A. Jain, A Wavelength Diver- 
sity Technique for Smoothing of Speckle, Ph.D. Thesis, 
California Institute of Technology, Pasadena, Califor- 
nia, 1975], [N. George and A. Jain, “Space and Wave- 
length Dependence of Speckle Intensity,” Appl Phys., 4, 
201-212, 1974], and [George, et al., “Experiments on 
the Space and Wavelength Dependence of Speckle,” 
Appl. Phys., 1, 157, 1975]. The surface heights obtained 
by varying the wavelength can be determined quite 
accurately to the order of a wavelength, the accuracy 
improving with greater heights. 

While the wavelengths of the ocean waves can be 
determined by inspection of the radar imagery or by 
obtaining its Fourier transform [W. E. Brown, Jr., C. 
Elachi, and T. W. Thompson, “Radar Imaging of 
Ocean Surface Patterns,” J. Geophys. Res., Vol. 81, No. 
15, pp. 2657 (May 1976); C. Elachi, “Two Dimensional 
Wavelength Spectrum of Swells Across the North At- 
lantic,” J. Geophys. Res., (in press)], it is also of consider- 
able interest to measure the heights of the ocean waves. 
A method of processing the images obtained by syn- 
thetic aperture radar to determine the heights of ocean 
waves at localized areas of interest in the images, and 
experimental results demonstrating this technique are 
presented. 

SUMMARY OF THE INVENTION 

In accordance with the invention, significant height 
information is obtained (of peaks in a rough terrain or 
heights of ocean waves) from a synthetic-aperture radar 
signal, S(0, by means for compressing the radar signal 
for intensity images h(x), l 2 (x) . . . l n (x) over different 
frequency spacings, A/*), of the available chirp or Dop- 
pler bandwidths, means for cross-correlating these in- 
tensity images over a selected area and means for deter- 
mining the value of a given (e.g., zero) component of 
this cross-correlation as a function of the spacing A/^of 
the center frequencies, A/o, for the cross-correlated 
images as a direct indication of height. In a digital pro- 
cessor, the data compression means is comprised of a 
plurality of match filters of predetermined width, A/„ 
each followed by a squarer to provide the intensity 
images I|(x), l 2 (x) . . . I n (x) which are cross-correlated 
over selected areas in the x dimension of the image to 
obtain the values /Ii(x)I*(x)dx for k selected from the 
sequence 1, 2, ... n. The cross-correlation value for k = 

1 yields a normalizing factor that may be used in subse- 
quent cross-correlations for k = 2, 3 . . . n. In an optical 
processor the data compression means is comprised of a 
data film on which the radar signal has been recorded, 
and a tilted plane correlator to produce on a map film 
the intensity images Ij(x), l 2 (x) . . . I„(x) for different 
frequency spacings, A/q. The cross-correlating means is 
then comprised of optical means for superimposing a 
selected area of an intensity image on a map film at one 
plane onto a corresponding area selected in the intensity 
images to be cross-correlated on a duplicate map film in 
another plane, and illuminating the superimposed areas 
with coherent light, and means for focusing the light 
transmitted through the superimposed and selected 
areas onto an intensity detector on a third plane. In 
either case, the apparatus thus includes means for com- 
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pressing synthetic-aperture radar data over different 
sections of the available chirp or Doppler bandwidths 
to obtain intensity images, and means for cross-correlat- 
ing these images to obtain height data. This height data 
for a predetermined range of frequency spacing, A/o, is 5 
plotted to calibrate the system. Once calibrated, only 
two intensity images need be cross-correlated at a pre- 
determined spacing to determine height directly. 


BRIEF DESCRIPTION OF THE DRAWINGS 


10 


FIGS. 1(a) and 1(b) show a system diagram for a 
height processor in which a synthetic-aperture radar 
signal S(0 is compressed by matched filters f\(t,x),f 2 (t,x) 

. . . f n (t,x), corresponding to a processing bandwidth of 
A/ r and at chirp or Doppler processing frequencies f^ x , 15 
fa 2 , . . . fan, respectively, the square of these outputs 
cross-correlated by the operation x, and the slope of the 
resulting curve measured to provide the height data. 

FIGS. 2(a) and 2(b) show radar imaging geometry 
useful in understanding the theory of the present inven- 20 
tion. 

FIG. 3 illustrates the configuration of a tilted plane 
optical processor, for compressing into a map film the 
radar signal from a data film, the bandwidth of the chirp 
or Doppler spectrum being controlled by the size of 25 
aperture A of a frequency plane filter F, and the center 
frequency being controlled by the position of aperture 
A. 

FIG. 4 illustrates the configuration of an optical 
cross-correlator in which the pupil P limits the area of a 30 
map film Ij at plane I illuminated by a laser L and colli- 
mator C, and lenses L n and L 22 of focal lengths 7(30.5 
cm each) image I lf onto plane II. The lens L 3 of focal 
length f Fourier transforms the product I t L 2 onto plane 
III, where the total power is detected by detector D. 35 

FIG. 5 illustrates cross-correlation data curves, 
R £ (A/o), for the ocean wave heights of 0.9, 1.1, 1.8, and 
4 m as a function of frequency separation for the chirp 
center frequencies used in the processing using a chirp 
bandwidth equal to 2 MHz. 40 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

Although reference will hereafter be made to the 
ocean surface in determining its roughness (wave 45 
height), it will be apparent that the same techniques 
may be employed to determine the roughness (peak 
heights) of land surface. As noted hereinbefore, the use 
of synthetic aperture imaging radars in oceanography 
has recently attracted considerable interest (W. E. 50 
Brown, Jr., C. Elachi and T. W. Thompson, Oceano- 
graphic Observation with Imaging Radar, URSI Fall 
Meeting, Boulder, Colorado, 1973; L. I. Moskowitz, 
“The Feasibility of Ocean Current Mapping vis. Syn- 
thetic Aperture Radar Methods,” Proc. Fall Conv. 55 
Amer. Soc. Photogrammetry, 760-771, Oct. 1973; T. R. 
Larson and J. W. Wright, “Imaging Ocean Current 
Gradients with Synthetic Aperture Radar,” URSI Fall 
Meeting, Boulder, Colorado, 1974). The radar generates 
a two-dimensional high resolution image of the ocean 60 
surface. The relation between the image brightness 
distribution and the properties of the surface is still not 
well understood; however three models (or any combi- 
nation of them) seem to be the most popular: 1. The 
image represents the modulation of the small scale sur- 65 
face roughness. This modulation could result from the 
presence of a swell, internal wave, local wind, slicks, . . 

. 2. The image represents the modulation of the local 


surface slope due to the presence of a large wave or ship 
wake. 3. The image is a result of the modulation of the 
water velocity by a wave or a current. The last model, 
if correct, is inherent to synthetic aperture imaging 
radars, in contrast to real aperture imaging radars, be- 
cause of the importance of the Doppler information in 
the generation of the high resolution image (A. W. 
Rihaczek, Principles of High Resolution Radar, McGraw 
Hill, N.Y., 1969; J. W. Goodman, Introduction to Fourier 
Optics, McGraw Hill, San Francisco, 1968; R. O. 
Harger, Synthetic Aperture Radar Systems, Academic 
Press, N.Y., 1970). Even though the formation of the 
image is not yet well understood, radar imagery does 
show a large variety of ocean surface phenomena 
which, in most cases, have been verified with ground 
truth measurements: large gravity waves, internal 
waves, currents, ship wakes, island wind shadows, sur- 
face slicks, weather fronts and others. It should be men- 
tioned that the uniqueness of the synthetic aperture 
imaging radar is that it provides high resolution imag- 
ery (tenths of meters) regardless of weather conditions, 
cloud coverage, sun illumination and platform altitude 
(even from orbiting spacecrafts). 

The radar image of long ocean waves provides two- 
dimensional patterns and therefore the wavelength of 
these waves. But an important parameter of these waves 
is their height. The system of the present invention 
shown schematically in the block diagram of FIG. 1 can 
be used to derive the wave height from the radar image 
by measuring the rate of change of the speckle pattern 
in the synthetic aperture radar image as a function of the 
radar frequency. 

A somewhat similar technique has also been used by 
D. Weisman, “Two Frequency Radar Interferometry 
Applied to the Measurement of Ocean Wave Height,” 
IEEE Trans, on Ant. and Prop., AP-21, 649-656 (1973) 
to derive ocean wave height along a flight line using a 
two-frequency correlation radar technique. 

In the art of radar imagery, a speckle pattern is de- 
fined as the pattern which results from the interferomet- 
ric addition of the echoes generated when a coherent 
electromagnetic wave is scattered from a rough surface. 
This pattern is familiar to people who observe surfaces 
using coherent laser light. (This frequency diversity 
technique has been extensively used in the field of co- 
herent optics to determine surface roughness.) If the 
surface is smooth there is no speckling. If the surface is 
slightly rough speckling appears. Then, a change in the 
electromagnetic wave frequency results in a change of 
the speckle pattern due to the change in the phase rela- 
tionships for the electromagnetic field scattered across 
the rough surface, the magnitude of this change increas- 
ing with the height variation on the surface. A detailed 
statistical analysis gives the resulting value for the cross 
correlation of the radar images of such a surface at the 
two different frequencies separated by Av, as 


{/ (v) I (v + Av)} - {/ (v)} 2 L. h o ^ f 

where { } denotes the average value, and I(v) is the 
intensity distribution in the image formed at frequency 
v, c the velocity of light, h 0 the rms height variation on 
the surface and /a factor depending upon the process 
used to obtain the imagery and the cross correlation, 
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i.e., film nonlinearities, changes from signal film to 
image film, image negative to image positive, etc. 

To implement the principle of above equation with 
side-looking synthetic aperture imaging radar 10 shown 
in FIGS. 2 a and 2b, the radar may have a chirped signal 
with frequency ranging from 1200 MHz to 1210 MHz. 
The echo is recorded on a data film S which may then 
be optically processed to generate map films as will be 
described with reference to FIG. 3 for use in an optical 
correlator shown in FIG. 4. Alternatively, the echo 
signal S(f) may be processed digitally to obtain wave 
height data as shown in FIG. 1. 

The principle involved in the ocean wave height 
measurement using side-looking (or forward-looking) 
synthetic-aperture imaging radar is that the resultant 
intensity for the electric field scattered from a rough 
surface at some arbitrary point in space changes as the 
wavelength or angle of illumination for the electromag- 
netic radiation varies. This variation is dependent on the 
height difference between the crest and trough of the 
wave, and a measurement of the rate of this variation 
provides an indication of the ocean wave height (or 
height of peaks in rough terrain). 

The essential elements in the processing of the syn- 
thetic-aperture radar data for obtaining significant wave 
height information are shown in FIG. 1(a). They in- 
volve compressing the radar data signal S(f) for ocean 
images over different sections of the available chirp or 
Doppler bandwidths/^, - • • /o«, cross-correlating 
these images and plotting the zero component of this 
cross-correlation value as a function of the spacings of 
the center frequencies for the bandwidths. A measure of 
the value of this curve for a given frequency spacing, 
A/„ would provide a direct indication of the ocean 
wave heights. This cross-correlation function, which is 
the output of the radar wave height processor shown in 
FIG. 1(a) is equal to the square of the product of the 
characteristic function of the ocean wave heights and 
the pupil function used in the chirp spectrum aperture 
for processing the radar signal. The intermediate calcu- 
lation of the image produced by the radar system, will 
take into account the effect of the orbital velocity of the 
waves when the optical system of FIG. 3 is used with 
the optical correlator of FIG. 4. A detailed treatment of 
this effect will be discussed hereinafter. Thus, a knowl- 
edge of the chirp spectrum pupil function of the radar, 
and an experimental determination of the cross-correla- 
tion curve by the wave height processor, leads directly 
to a measurement of the ocean wave height distribution. 
It has been found that the measured cross-correlation 
curve is equal to T XA/k)) exp (—0.1 where cr H 

is the significant wave height, A/^the frequency separa- 
tion for the cross-correlation value, and T/A/o) a num- 
ber depending upon the chirp spectrum aperture func- 
tion I%for the radar processor. 

Referring now to the processor shown schematically 
in a functional block diagram in FIG. 1(a), a number of 
functional blocks 1, 2 ... N receive the chirp return 
signal S(l) and split it into n frequency channels, each of 
a predetermined width, A /„ of about 2 MHz. In prac- 
tice, digital bandwidth filters may be used to provide 
the n frequency channels for electronic data processing, 
but in the preferred embodiment to be described with 
reference to FIGS. 3 and 4, optical filtering techniques 
are used. Outputs from the functional blocks 1 through 
N are individually squared in circuits 11, 12 . . . NN to 
change the electric field signals to intensity signals I t (x), 
l 2 (x) . . . l n (x). These signals are cross correlated in pairs 
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by a circuit 20 and integrated. A slope detector 21, 
calibrated in heights, detects the drop-off (width) of the 
resulting integrated signal to provide as a system output 
a direct measurement in height. This technique for 
5 height determination is based on the fact that the syn- 
thetic aperture radar is an imaging instrument, where 
the resolution is determined by the chirp and doppler 
bandwidths used for the radar system. Thus, any point 
on the image contains contributions from the field scat- 
10 tered by an area of the terrain determined by the resolu- 
tion cell size of the radar. By narrowing the effective 
doppler or chirp bandwidth utilized by the instrument, 
the area of the ocean surface contributing to the inten- 
sity at a given image point may be widened to any de- 
15 sired magnitude. The total path length that the radar 
wave has to travel to the image is different for the 
trough of the wave than for the crest, and the resultant 
phase difference for the contribution from these two 
regions of the wave is dependent on the frequency or 
20 angle of illumination of the wave. As this frequency or 
angle of illumination is varied, the fluctuation of the 
intensity at the image point, which receives contribu- 
tions from both wave regions, is related to the differ- 
ence in heights between the crest and trough of the 
25 wave, modified by the effect due to the resolution cell 
size of the radar. By measuring the effective rate of this 
fluctuation, one can infer the heights (of the ocean 
waves or of peaks in a rough terrain). 

Before processing with a more detailed description of 
30 wave (or peak) height determination, the optical system 
shown in FIGS. 3 and 4 will be described. It has been 
used with an L-band synthetic aperture imaging radar 
operating at 25 cm wavelength, with a phased array 
antenna 75 cm long, 25 cm wide, and 5 cm thick to 
35 obtain a range beamwidth of 90% centered 45° off verti- 
cal and an azimuth beamwidth of 18% centered on zero 
doppler. The transmitted peak power is 7 kW, the pulse 
repetition frequency 1 kHz, the radar platform velocity 
250 m/sec, the pulsewidth 1.25 jusec, and the total chirp 
40 spectrum bandwidth 10 mHz. The data are recorded 
optically on a film, shown as a data film 22 in FIG. 3 
using a dual-channel recorder, such as the Goodyear 
102, dual-channel recorder, and processed optically to 
convert the signal film into imagery with a resolution of 
45 20 m range (45°) and 10 m along track. Once the data 
film is developed, it is employed in a tilted plane corre- 
lator (TPC) shown schematically in FIG. 3. The data 
film 22 is illuminated by the output of a HeNe laser (not 
shown) after it has passed through a spatial filter-col- 
50 limating lens system. A Fourier transform lens L lf a 
frequency plane aperture filter F, and an imaging lens 
^convert the range data into the range image at a plane 
on a map film 23. A cylindrical azimuth telescope L 3 has 
the function of bringing the azimuth data at plane I to an 
55 alongtrack image focus on the map film 23. By control- 
ling the position of a data plane aperture 24 and the 
spatial position of the frequency plane aperture filter F, 
the chirp frequency and the center bandwidth can be 
controlled. By setting the value of the frequency plane 
60 aperture at a given range bandwidth, the same strip of 
signal film can be processed for different center fre- 
quencies to provide the required images by positioning 
the location of the frequency plane aperture at the cor- 
responding spatial coordinate in the range dimension. In 
65 that manner, the radar return signal recorded on the 
data film from a synthetic aperture radar in an aircraft 
30 (FIG. 2a) is processed by the TCP for different chirp 
bandwidths . ./on of predetermined width (e.g., 
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2 MHz) to provide on the map film 23 a record of the 
split radar signal equivalent to that of the different chan- 
nels 1, 2 ... n in FIG. 1(a). The map film itself functions 
as a squarer. 

The optical setup used to cross correlate pairs of 5 
images of the same area, but obtained for different chirp 
spectrum center frequencies, is shown in FIG. 4. An 
argon-ion laser 31, combined with a spatial filter colli- 
mator system, C, is used to illuminate the film image 1Q 
I, (*, y, fa) processed for the particular chirp spectrum 
frequency fa and the bandwidth f r The two Fourier 
transform lenses L u and L 22 , each of focal length/ are 
used to provide unity magnification of an image at plane 
I onto plane II, where the transparency of transmissiv- 15 
ity I 2 (x + A x,y + Ay, fa + A/o) is situated, which is the 
image of the same area as Ii(x, y, ff) but processed at the 
chirp spectrum frequency fa + A/oand shifted by (Ax, 

Ay) in the X- and y-directions for the optical system. 
The distance between L n and L 22 is 2/ and between 20 
plane I and L n and plane II and L 22 is / so that no phase 
effects are introduced by the imaging process. The 
spatial distribution of the electric field at the exit of 
plane II is proportional to Ii(x, y, f^l 2 (x + Ax, y + Ay, 25 
fa + A/o), and this is Fourier transformed by a lens L 33 
of focal length / onto plane III. The power incident on 
plane III is measured by a photodetector 32 whose 
aperture is large enough to collect all the Fourier com- 
ponents involved. Using Parseval’s theorem, and assum- 30 
ing a unit amplitude illumination at I lf the power mea- 
sured by the detector is 


4*0 .Vo) = If »K*o JV * >04* y)dxdy . 
- 00 


( 1 ) 


W^Af^Ax.Ay) = 


35 

jfjf \Ii(x,y,f K Mx+Ax, y+Ay, f^Af^dxdy 


where A/ is the bandwidth at which and I 2 were pro- 
cessed, (Ax, Ay) are relative shifts of Ij with respect to 
I 2 , and fa fa 4- A/o are the two processing frequencies 
involved. For the purposes of wave height measure- 
ment, Ax and Ay would be equated to zero and the 
function Wvr(A/ro, 0, 0) plotted as a function of A/o. In 45 
considering Ij and I 2 to be the two images of the ocean 
wave, we have assumed that the film nonlinearities of 
the image film can be ignored and have factored out the 
gamma of the film. The result is displayed using a CRT 
33 or recorder 34, and a slope detector 35 provides a 50 
direct indication of ocean wave height after the system 
has been calibrated. 

As noted hereinbefore, as the frequency or angle of 
illumination is varied, the fluctuation of the intensity of 
the image point, which receives contributions from 55 
both wave regions is related to the difference in heights 
between the crest and trough of the wave, modified by 
the effect due to the resolution. A detailed discussion 
will now follow on the theory by which heights, of ^ 
ocean waves or peaks of rough terrain, can be inferred 
from measurements of the effective rate of this fluctua- 
tion. 

In the linear system idealization for the synthetic 
aperture radar performance [Radar Handbook, ed. Mer- 65 
rill I. Skolnik, Chap. 23 (by L. J. Cutrona), New York, 
1970] the electric field e(xo y 0 ) in the image of the radar 
output is 


wherein y <>* x, y) is the generalized ambiguity func- 

tion as defined by Cutrona [ibid.] for the radar and g(x, 
y) is the terrain reflectivity. 

In the calculation for the ambiguity function, how- 
ever, it is assumed that the terrain reflectivity does not 
depend upon the wavelength of illumination or the look 
angle of the radar. In the case of an uneven surface, this 
assumption is not applicable, and it is necessary to con- 
sider the effects of the wavelength and angular depen- 
dence of the radiation scattered by the surface. In order 
to take this phenomenon into account, calculate the 
image field for each wavelength component of the chirp 
spectrum and position of the radar platform and inte- 
grate this electric field over the total synthetic antenna 
length and the chirp spectrum bandwidth available. 

The impulse response function for the radar imaging 
system for a given value of the antenna position, x A , and 
the chirp spectrum frequency,/, may be expressed as a 
convolution of the generalized ambiguity function with 
the Fourier transform of the pupil function for this 
range and Doppler aperture. This aperture is the Diarc 
delta function 8(f dd — f d )8(f rr — /), where the phase 
factor due to he time delay, f d , is equal to xJ'a/Z , fa 
being the frequency of operation of the radar and Z the 
height of the radar platform from the terrain. While the 
Doppler shift of the return signal is fivr/c, for v r the 
aircraft velocity and c — 3 X 10 8 m/sec, it is the quan- 
tity / d that is significant in the matched filter processing 
of the radar correlator. The Fourier transform of this 
delta function is defined by ' 


u{*> yo> - 

ff 8 (fo - - />*p £ 


( 2 ) 


—fly 


(fid* + M 


~\d}^f rr 


the symbol.^ being used to specify the particular nature 
of the chirp processing operation. 

The resultant impulse response of the radar system 
for these particular frequencies (fo /) is given by the 
convolution of the impulse response U(x, y Q ) with the 
generalized ambiguity function \jf(x, y ). Thus, for the 
complex terrain reflectivity g(x, y) t the electric field 
recorded in the radar image is the convolution of g(x, y) 
with the convolution of x/j(x, y) with U(x, y 0 ). Further, 
for the Doppler spectrum aperture function given by 
rff d — fdo)> which has an effective width of Af d and is 
centered at fa, and a chirp spectrum aperture function 
Tff r — /o) of width A/ centered at the radar frequency 
fay the electric field in the radar image is the superposi- 
tion of fields obtained for each Doppler and chirp fre- 
quency component weighted by the functions Tff d — 
fa), r ffr — fa). This yields the equation 


4*o yd) = 


- ri 


r -jv (3) 

g{x, y)*\fi(x, J.)*exp {fix + M \\ 

X r Jf,,- 


where * denotes the convolution operation. Eq. (3) 
describes the field in the image plane of the terrain 
reflectivity g(x, y) obtained by a synthetic aperture 
radar defined by the generalized ambiguity function 
i|f(x, y) and the processing Doppler and chirp spectrum 
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apertures TJf d — /o), Tff r — f K ,). Since the azimuth 
aperture for the synthetic aperture radar processor is 
generally symmetric, it is safe to consider the value of 
fdo to be zero in our further evaluation of Eq. (3). 

From the geometry of FIG. 2 a the terrain reflectivity 5 
function gp( x, y) for a particular area surrounding the 
point P located at (0, y p ) is: 


10 

-continued 


SM y) = R ( x > y) exp [-~/,cos Qpxyfri** y)] exp^-“- 


(4) 
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X fXx sin(0 A - 20^ + y sin(a^ - 26^) 
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R(x, y) is the amplitude reflectivity of the terrain due to 
absorption and dielectric constant variations and h(x, y) 
is the profile of the terrain resting on the local slopes of 
angles ( $ P w Bp^) shown in FIG. 2b for the radar plat- 
form and terrain of FIG. 2a which vary slowly com- 20 
pared to the variation of h(x, y). The slope angles 0^ 
and ^are a function of both the slope angles ( 0 R „ Bpy ) 
of the side-looking radar and the tilt ( t ) of the terrain 
which, for an ocean wave, is the angle between a tan- 
gent to the ocean wave and a horizontal X-Y plane. 25 
FIG. 2b illustrates the geometry of the angle Bp^ mea- 
sured in the Y-Z plane and FIG. 2c illustrates the geom- 
etry of the angle Bp^ measured in the X-Z plane. Note 
that in FIG. 2 a there is a slight tilt ($& 6#) of the terrain 
with respect to the X and Y axes of the horizontal X-Y 30 
plane, to illustrate a typical case in which the roughness 
of the terrain is the major factor in the tilt a ngles Bp^ 
Bpa. The distance of point P from the radair, VZ 2 +y p \ 
is considered to be large, so that the radar wave incident 
on the area surrounding P is essentially plane and is at 35 
angles (B^ Bpy) with respect to the normal to the hori- 
zontal plane. The function exp[(/'47r/c)/cos0 y)] is 
the phase modulation in the reflected wave due to the 
small-scale roughness h(x, y). The factor 


( air 

— 


fix sin(0 ft - 20ft*) + y sin(0^ - 20^)] 


) 


is the phase modulation in the reflected wave from the 
tilt variations of the ocean wave and the angles of inci- 
dence Bpx and 0^of the radar wave and is determined 
from equating the angle of incidence to the angle of 
reflection from the surface. The angle B Px is generally 50 
small, and for a typical spacecraft radar it is less than 0. 1 
radian, (fip^ Bp^) are the angles of tilt of an ocean wave 
when g p (x, y) is an ocean surface and are generally lim- 
ited to be less than 0.143. Thus, in further computations 
involving gj(x, y) it is reasonable to consider sin(0/* — 55 
20/* x) as approximately equal to (sin0/* — SU 120 P**). 
Since sin0/* is equal to *yz, /osin0/>* is the Doppler 
frequency f d , and 0/^is generally much larger than B^y 
for side-looking radars. 

To evaluate the radar performance for a given reflec- 60 
tivity, substitute Eq. (4) into Eq. (3). For h(x> y) = Bp^ 

= Bpty =s 0, which corresponds to an even terrain of 
reflectivity R(x, y), the image electric field is found to 

be ■■ ■ 

65 


^ 2x 0 Pq sin Opy ^ 


15 Eq. (5) neglects the variation of x sin 0/*in comparison 
to y sin 0/^ which is a much larger term, in the function 
yXy)* For infinite processing bandwidths, given by the 
widths of T X/p> T/fJ), the functions TX*), y/y) reduce 
to delta functions, and Eq. (5) reduces to Eq. (1). 

On the other hand, for a radar with an infinitely high 
resolution, when xf/(x, y) would be a delta function, the 
image electric field would be the convolution of the 
reflectivity with the impulse response functions T rf and 
T„ and the image resolution would be determined pri- 
marily by the bandwidths utilized by the radar proces- 
sor. In general, the resolution of the radar is indeed 
limited by the data handling capabilities of the proces- 
sor, and so it is safe to assume, for further analysis, that 
the generalized ambiguity-function is the delta function. 

Consider the case of an uneven surface, such as ocean 
waves, and assume that the function R(x, y) is unity in 
Eq. (4). We describe h(x, y) as a series of steps whose 
perimeters are defined by the functions S^x — x? y — 
y) centered at the various points (x> yj) 9 and the profile 
of the surface within this (ij) th step is defined by the 
function h ip The effective length of the S, y function is 
given by in the ^-direction and lj in the ^-direction. 
Thus, the surface profile h(x, y) can be written as 

2 hijSij(x - x,y - y). 

U 

For nonoverlapping S iJf the exponential exp [(—/477/c) 
cos Bpxyhix, y)] can be rewritten and the resulting sur- 
face reflectivity can be written in the following form: 

f \ (6) 

gfay) = exp { sin (0/* - 20 Pu) + y sin(0^ - 2dpf y )] \ 

x{l + ZSifoc -x,y- J^exp ^-^-/ r cos 0/*, AyJ- 1^ 

The first term in Eq. (6) represents the field reflected if 
the terrain is perfectly flat with only tilt variations. The 
term identified by the summation consists of two terms, 
the first being the component of the reflected field due 
to the surface roughness, and the second being the com- 
ponent that must be subtracted from the field for the flat 
surface case to provide the reflected field for the 
smooth component of the terrain. 

To obtain the electric field in the image of the field 
given by Eq. (6), substitute Eq. (6) in Eq. (3), neglect the 
effect of the variation of //compared to f d in the term 
cos0^) and assume that the scattering cross-section of 
( 5 ) the terrain, determined by the Fourier transforms of 
Sifpc — x ir y — yj) does not vary significantly with the 
Doppler chirp spectrum frequency variations over the 
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bands used in the radar processor, as compared to the 
variations of the exponential phase terms containing this 
(fd, ft dependence. The resultant image electric field is 


12 


-continued 
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where fly is the scattering cross-section factor and is 
equal to 

sin 2<W sin(^ - 29 ^ ] 

for Sy defined by 


**«>- ff 


Sf/x,y) exp [i2ir(x! + ym)] dxdy. 


30 


In the approximation, where S t pc, y) are given by the 
step functions [rect(x/lj)], the cross-section factor Sy is 
equal to 

35 



lj sin 2 Opt* 



The functions y /2x/c), y Jy sin(6py — 2 dp#)] are the ^ 
same as those defined in Eq. (5) and, for the chirp and 
Doppler spectrum pupil functions, approximated by 


rjfj = reel (-^-)and TfQ = reel 45 

these functions become 

A/i sine ^ J A/, sine j^T" y sin (6^ - 26 hy ) J. 

The product %*y r on the left-hand side of the convolu- 
tion sign in Eq. (7) corresponds to the impulse response 55 
of the radar, and the resolution cell size is determined to 
be 2c/[bf r sm(Qpy — 26^)] in the range dimension and 
c/A/rf in the azimuth dimension, this being equal to the 
size of the radar antenna if the full Doppler aperture is 
used. The function F ,/jc, y; x ¥ yf h 0 ) describes the effect 60 
on the image field due to the phase variations in the 
scattered field arising from the unevenness of the terrain 
and is given by 

Fifx> y; Xf y; h ( ) = 65 

S . */sin 20ft* — A# cos % 1 

sin (Spy — 26^) J~ 


where cos Opyo is equal to the value of cos 0py X at the 
Doppler frequency f d equal to zero. 

The scattering factor ft is dependent on the tilts of the 
ocean waves (Op^Op^) and also the distance away from 
nadir as determined by 0^and the size of the small scale 
roughness (1„ l y ). A modulation in the wave tilt or 
roughness size would contribute to a modulation of the 
intensity in the image. A detailed analysis and experi- 
mental evaluation of the effect of tilt and roughness of 
the terrain on the image modulations for the radar sys- 
tem have been reported by Schaber, et al. [G. G. Scha- 
ber, G. L. Berlin, and W. E. Brown, Jr., Geol. Soc. 
Amer. Bull , March 1975] for the case of desert images 
and by Elachi and Brown [C. Elachi and W. E. Brown, 
Jr., IEEE Trans, on Ant. and Prop., in press] for the 
case of ocean wave imaging. 

While in principle, the synthetic antenna length is 
determined by the area illuminated by the radar, the 
actual length is modified by the finite angle over which 
the radiation is backscattered by the small-scale rough- 
ness. Since this angle is determined by the spatial 
lengths of these scatterers, a modulation in the size of 
the scatterers would also result in a modulation in the 
azimuth resolution of the radar if the scatterer sizes 
were larger than the size of the radar antenna. 

Note that Eq. (7) does not include the effect of mo- 
tions of the surface being imaged. Since f d is the time 
variable in Eq. (3), the effect of surface motions can be 
included by making g(x, y) a function of f d as well. A 
uniform motion of g(x, y) in the y-direction with a ve- 
locity v g can be included by substituting the function g(x, 
y - Zfd) for gix t y) in Eq. (3). £ is equal to Zv g /ftv r \ 
where v r is the velocity of the radar platform. A result- 
ing evaluation provides no change in the image except 
for a uniform shift proportional to £ in the azimuth 
direction. A uniform motion in the x-direction is equiva- 
lent to a change in the equivalent velocity of the aircraft 
and thus a change in the focusing parameter for radar 
imagery. 

At any one point on the ocean surface, the height 
varies according to the orbital velocity, and the cross- 
section of this point according to the ocean wave 
length. A patch of constant cross-section moves in the 
direction of the ocean wave propagation at the ocean 
wave velocity. Since the radar maps the time delay of 
the backscattered signal from a patch of constant cross- 
section, the ocean wave appears to it to be moving in a 
plane parallel to the radar platform motion and the 
orbital velocity effects are not important. However, an 
object of constant cross-section, such as a ship, resting 
at a particular point on the ocean, and moving vertically 
at the orbital velocity, will have a modified image and 
its image will generally be degraded in the azimuth 
direction. It should be pointed out that the synthetic 
aperture radar is sensitive to the variation of time delays 
of the signal scattered by a target, and the actual effect 
of the instantaneous shift in the frequency of the back- 
scattered signal due to the orbital motion is included in 
this time delay variation which depends on the variation 
of position of the scattering patch. In such a consider- 
ation, the variation of the velocity and radar cross-sec- 
tion of this patch must also be taken into account, and 
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approximating these quantities by a constant number 
during the integration time of the radar processor has in 
the past led to considerable confusion in the interpreta- 
tion of the synthetic aperture radar ocean wave imag- 
ery. A dependence of the ocean wave tilt angle 6^ on 5 
time would lead to a change in the resolution cell size of 
the radar, the result being a smaller resolution cell size 
if the tilt angle is increasing and a larger cell size if it is 
decreasing. This effect, however, would be very small 
for the normally occurring rates of change of the tilt, \q 
The random motion of the small-scale roughness would 
have negligible overall effect on the radar image since 
its main contribution is in the backscatter modulation 
function fi ip a treatment of this negligible effect has been 
considered by Elachi and Evans [C. Elachi and D. E. 15 
Evans, IEEE Trans, on Ant and Prop,, in press]. 

The ability to determine height of ocean waves, or 
peaks of a rough terrain, according to the invention 
derives from Eq. (7). For a particular point Pi in the 
area of interest, the image electric field may be rewrit- 2 o 
ten as 
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e p i (xp l > y Pv /*)) — 

exp {^-f*>y Py ^ 

Hit 


( 8 ) 


25 


exp 


£ — Op, sin ‘IQptyj 4- *,sin 2 Op*) J 


corresponds to the phase variation due to the variations 
in the tilt angles of the terrain. If the resolution cell is 
small enough so that the tilt of the waves does not 


change appreciably across its length, the decorrelation 
length due to the tUt in the range and azimuth directions 
is determined by A//(2 sinfl/^) X (sinflp/cosflp,) and 
A//(4 sin#/**), respectively. The composite decorrela- 
tion length would depend upon the sum of the slopes of 
the individual decorrelation curves for the resolution 
cell size and the tilts approximated by 


4 f r 


+ | fry exp ^^f l0 (y i sm(e Py - 2 9^ - *,sin 20 m pJ 

X £exp cos 0^ h^-X 3 ° 

where the finite size of yfy), yj(x) limits the number of 
scatterers that contribute to the field at an image point. 

In this case, this number is given by the summation over 35 
the N scatterers; a is the uniform background level term 
and is the first term on the right-hand side of Eq. (7); fi 
is the scattering factor; and (0^ dp^) are the tilt angles 
at points (x, y) which have been assumed to be slowly 
varying compared to the small-scale roughness h(x, y), 40 
Assuming that the scatterers are densely packed [L. I. 
Goldfischer, J, Opt Soc, Am,, 55, 247-253, 1965], N is 
equal to the resolution cell area divided by the area 
occupied by a scatterer, and this quantity can be ap- 
proximated by 2c 2 -/(l l l j Af r hf d sin Opy) (neglecting the 45 
effect on the resolution cell size due to the small tilt 
angles, the scatterer size modulations, and the heights 

h). 

The small-scale roughness is usually less than the 
wavelength of the radar and the variation of ex- 50 
p[(/4n*/c)/Jo cos Opycfiij] with respect to fa is negligibly 
small over the chirp spectrum bandwidth available. 
Also, since 0^ is a small number, it is permissible to 
approximate sin(0 py - 20^) by (sin 0^- cos 0p,sin 26^); 
The term exp[(i27r/cy^/sin Op) represents the vana- 55 
tion of the phase of the radar wave incident at angle Opy 
across the resolution cell size of the radar, and with 
respect to change in radar frequency, occurs in an inter- 
val equal to the chirp bandwidth of the radar processor. 
The factor 60 


1 + 2(cos0p/sin0p y ) sin0^ + (4 sin dp^/Afy ' 

The resolution cell size of the radar may be large 
enough to cover more than half an ocean wave, in 
which case the tilt angles (Op^ Opt) cannot be assumed 
to be constant within a resolution cell of the radar, but 
it is possible to rewrite (y ; cos fl^sinfl^. -1- x,sin20^ as 
2 H $ corresponding to the height of the ocean wave. 
Due to the cos Opy term, the dependence of the resolu- 
tion cell size on sin 0/^and the tilt angle variations, the 
Hy as defined is not strictly the ocean wave height pa- 
rameter, but for the purposes of the statistical analysis, it 
may be considered to be equivalent. Then calculate the 
cross-correlation for the image intensities processed for 
the two frequencies fa and (fa + A/q) and obtain the 
following result 

X \ (9a) 

+ 4/*,o) fa 4 + 4Afa 2 /? 2 + 2/J 4 iV(l + 27V)] + 

<t h ) ** )+ 1 ]!W*> + 

0W - ** )+ 1 ]r?(AW 

The brackets < > denote the mean value operation. /3 
is the average value of /^ within the resolution cell and 
is the backscattering cross-section of the area of interest. 
N; the number of scatterers in a resolution cell, is ap- 
proximately equal to 2c 2 /(AfAf d $in0 / > y I i /) t where If- are 
the sizes of the small-scale scatterers. cr^is the standard 
deviation of the ocean wave height function H ip and cr h 
is the standard deviation of the heights of the scattering 
waves. Tjtfand ij A are the characteristic functions for the 
probability density of the wave height H f y and the 
roughness h 0 and, for a Gaussian distribution, are equal 
to exp{ — \ [477 ( A/q/ c) <r H ] 2 } and exp{-|[47r(A/o/c) 
cr ft] 2 }, respectively. For the special case when the reso- 
lution cell size is smaller than an ocean wavelength, the 
characteristic function for. the wave height 0 H in Eq. 
(9a) may be replaced by the characteristic function for 
the wave slopes equal to 

Veti&AfrtAy cosfl^sin 6^) ^JlA/^AxO^) 

for Ax, Ay, the resolution cell sizes. For a Gaussian 
density function for the slopes, this function is equal to 


exp< 


1 n 


('cos« f> , \ -1: 



sin Opy 
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where <r etyf o-^are the standard deviations of the wave 
slopes in the y- and x-directions, respectively, 

a, which corresponds to the specular reflection term, 
is negligibly small compared to 0 for the non-nadir 
angles and may be neglected in Eq. (9a). Further, N is a 
very large number, and the terms of order N can be 
neglected in comparison to the terms containing N 2 . cr h 
is a small number compared to cr Hi and 'q h [2hf r0 /c)crf] 
may be approximated by unity for the values of A/o 
involved with the radar chirp. The quantity 

(iMW* + vS> -^-r - * 

is defined by R(A/o) as the normalized cross-correlation 
function for the radar image. This quantity is simply a 
product of the square of the radar processor pupil func- 
tion, r, 2 (A/o)> and the square of the characteristic func- 
tion of the wave heights or the characteristic function 
for the wave slopes. This may be written as 

/ \ (9b) 

*< Att-* 

for the radar resolution cell much larger than an ocean 
wavelength, and 

f w v W 

R (4/ri) “ rfey JT 2 - AV COS 9p/T api) J 

^«x( — 7° b xar $ptt (i/,o) 

for the radar resolution cell smaller than an ocean wave- 
length. Thus, an experimental determination of the 
cross-correlation function along with a prior knowl- 
edge of the pupil function in range for the radar proces- 
sor would lead to a determination of the characteristic 
function for the ocean wave heights from which the 
wave height spectra of the ocean waves may be deter- 
mined. When the radar resolution cell is larger than the 
ocean wavelength, Eq. (96) can be used to interpret the 
radar data to obtain information on the distribution of 
wave heights. When the radar resolution is smaller than 
an ocean wavelength, Eq. (9c) provides a distribution of 
wave slopes, and knowledge of the radar resolution cell 
size can be used to convert this information into wave 
height information. Also, note by examining the rate of 
change of R(A/o) as a function of cr H , that the sensitiv- 
ity of the measurement of cr H from R(A/ K) ) increases 
with o> The pupil fraction T/A/o) in Eq. (9a) has the 
effective width (A//A/q). For A/o larger than A//2, the 
cross-correlation function becomes very small, and for a 
measurement of the characteristic function t\ H , it is 
necessary to measure R(A/o) for A/o< A//2. While the 
theoretical form of the pupil function TXf r ) is a step 
function of bandwidth A /, practical radar systems do 
not realize this functional form because of the antenna 
pattern and nonlinearities in the data recording and 
processing operations. 

The result in Eq. (9a) could have been obtained di- 
rectly from simple physical considerations by adding up 
the phases of the electric field reflected by the terrain 
within a resolution cell area of the terrain. The pupil 
function factor TXA/)o) is a direct consequence of the 
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nonvertical angle of incidence of the radar radiation, 
and the characteristic function factor t} h is a conse- 
quence of the phase variations for reflection from a 
rough surface. 

In order to obtain an experimental evaluation of 
<h x (frQ0p x {f*> + A/o)> in Eq. (9a) for a given ocean 
area, the optical processor of FIG. 3 has been used to 
obtain the successive images I(jc, y, /o), I(x, y, /o + 
A/oi), • . . I(*, y fro + A/o rt ) at the frequencies/^ /<, + 
10 A/oi> . . ./o + A/o„ but for the fixed processing band- 
width of A f r The image l(x, y, /*,) in the map film 23 is 
placed in plane I of the optical cross-correlator in FIG. 
4, and using the time codes in these images as a refer- 
ence for accurate superposition of the image I(x, y, /o) 
5 on I(x, y, /o 4- A/**) at plane II, the quantity W Vr (A fa, 
0, 0) is measured for each successive l(x, y, /<, A/o*). 

In order to account for the effects of laser power fluctu- 
ations and background exposures on the film, a mea- 
2 q surement of W A / r (A /o*) is also made for large Ax, Ay, 
and this is subtracted from W^ r (A/o*, 0, 0). To normal- 
ize the effects of average intensity variations for differ- 
ent I(x, y, /o -f- A/o*) further, the difference 

25 IW V ,(4 fm, 0, 0) - A^ Ay)] 

is divided by the correlation number when there is no 
frequency shift, i.e., by W A/r (0, 0, 0) — W A/r (0, Ax, Ay). 
The resulting curve R £ (A/o)> which is defined by 

30 

^(A/^0,0) - WJW^bcAy) 

- »V r (0,Ax,A y) 

is plotted as a function of A/o and corresponds to the 
35 normalized value of | <I/> 1 (/m))I j p } (/>' / o + A/o)>| 2 , the 
area on the image correlated bemg controlled by the 
pupil P. The quantity R £ (A/o) measured by the cross- 
correlator is equivalent to R(A/o) defined in Equation 
(9a), the superscript E being used to differentiate be- 
40 tween the experimental and the theoretical functions 
involved. 

FIG. 5 is a plot of the function R^(A/o) from experi- 
mental data where the bandwidth used in processing the 
successive images was 2 MHz, the processing chirp 
45 frequency interval for each successive pair being 0.5 
MHz. Since the curves are a product of the characteris- 
tic function for the ocean wave heights and the radar 
pupil function, and since the radar pupil function falls 
50 quite drastically for the frequency separations greater 
than half the chirp bandwidth for the radar processor, it 
is necessary to confine the examination of the experi- 
mental curves for A/o < 1 MHz to obtain the wave 
height information. For a given ocean area, the curves 
55 were obtained by performing cross-correlations over 
different portions of that area, and the deviation of the 
values of R £ (A/o) for a given A/o was approximately 
±0.05. As expected, no effect was noticed for the varia- 
tion of R £ (A/o) for different positions along the y-axis 
60 since the radar aperture function is independent of the 
look angle. 

An examination of R £ (0.5 MHz) shows a variation in 
the wavelength decorrelation with respect to the ocean 
wave heights involved, and the points appear to fit the 
65 curve 0.65 exp(—0.066or w ), where cr H is the significant 
wave height. The coefficient 0.066 in the exponential is 
due to the decorrelation factor from phase changes 
from crest to trough of a wave, and is close to the theo- 


retical value of [167r(A/o/c)]. The factor of 0.65 is the 
value of the radar chirp spectrum function at 0.5 MHz. 

The dependence of the decorrelation function 
R £ (A/,o) on the chirp spectrum bandwidth used was 
verified by obtaining the curves R £ (A fo) for 6-, 4-, 2-, 5 
1-, and 0.5-MHz bandwidths, and extensive decorrela- 
tion was noted for the frequency separations A/o equal 
to half the bandwidth employed in the processing. For 
an infinite-bandwidth radar, the effect on the decorrela- 
tion should be nonexistent. By extrapolating the R £ (0.5 10 
MHz) for an infinite bandwidth and comparing this 
result for R £ (0.5 MHz) at 2 MHz, the value of 0.65 for 
the chirp aperture effect was verified. 

When the cross-correlation function, R £ (A/<>), is eval- 
uated for the frequency separation of A/o — 1 MHz for 15 
the different wave heights, the functional value ob- 
tained is 0.29 exp (0.12 <r H ). The doubling of the factor in 
the exponent is due to the doubling of A/o. The effect of 
the chirp spectrum aperture function is given by the 
factor 0.29 for this case. 20 

The chirp aperture function is the Fourier transform 
of the radar resolution in range. This resolution can be 
readily measured in the optical processor of FIG. 4 by 
plotting R £ (0) for different position separations of the 
two films being correlated, and an experimental exami- 25 
nation of this resolution function shows that the width 
of the resolution cell is inversely proportional to the 
chirp bandwidth used. 

In these results, attention has been focused on mea- 
suring the effect on the radar images due to different 30 
effective frequencies of the radar operation. A similar 
calculation and experimental procedure can be em- 
ployed to exploit the effect on the radar images of dif- 
ferent look angles for the radar, and this can be done by 
processing the radar data for different Doppler band- 35 
widths and cross-correlating these images as a function 
of the center Doppler frequency separations for these 
bandwidths. 

To summarize the operation of the digital processor 
shown in FIG. 1(a) in accordance with the system the- 40 
ory discussed in detail above, a synthetic-aperture radar 
return signal S(0 of bandwidth f r is separated into a 
plurality of frequency channels by functional blocks 
(match filters 1, 2 ... N, each channel of predetermined 
width, A f n over the entire available bandwidth as 45 
shown in FIG. 1(b) for a chirp or Doppler spectrum. 
The chirp spectrum is used for images to one side of the 
radar track, and the Doppler spectrum is used for im- 
ages along the radar track. The outputs of the match 
filters are electric field images which are squared by the 50 
squarers 11, 12 . . . NN to obtain intensity images I,(x), 
I 2 (jc) . . . I„(x). These images are degraded since only a 
part of the chirp or Doppler spectrum is used for each, . 
but height information is the objective, not a high reso- 
lution image. 55 

These intensity images are cross-correlated by the 
circuit 20, but first the image Ij(x) is auto-correlated 
(i.e., cross-correlated with itself) to obtain a normalizing 
factor to be applied to subsequent cross-correlations. 
Assuming a normalizing factor of unity for simplicity of 60 
description, the image Ij(x) is then cross-correlated 
with each of the images I 2 (x) through l n (x) over a se- 
lected area. These cross-correlation functions are calcu- 
lated and fixed (e.g., zero) components of these cross- 
correlation values are plotted as a function of the spac- 65 
ings of the center frequencies for the channels to obtain 
the characteristic functions of the ocean wave heights 
as shown in FIG. 5. These plotted curves calibrate the 
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processor for subsequent wave height measurements. 
Thus a measure of the value of these curves for a given 
frequency spacing, e.g. A/* = 0.5, provides a direct 
indication of the ocean wave heights as a function of the 
slope of the curves. Consequently, the output of the 
slope detector 21 is a direct indication of ocean wave 
height. In other words, after once calibrating the system 
from data plotted as in FIG. 5, the system will provide 
wave height data directly from the ratio of the cross- 
correlation function, R £ (A/o), to the frequency spacing, 
A/*. 

The functions of the digital processor shown in FIG. 
1(a) are implemented optically in an analogous manner. 
The optical processor shown in FIG. 3 performs the 
functions of the match filters 1, 2 . . . N of FIG. 1(a), and 
the map film, which records the electric field images 
thus derived from the radar signal recorded on a data 
film 22, inherently squares the electric field images thus 
performing the functions of the squarers 11, 12 . . . NN 
of FIG. 1(a). Two copies are made of one of the inten- 
sity images thus recorded, such as for channel f m . One 
is placed in plane I in the optical processor shown in 
FIG. 4 for cross-correlation with the other intensity 
images recorded on the other map film placed at plane 
II of FIG. 4. The pupil P selects the area of the intensity 
image in plane I to be cross-correlated. The cross-corre- 
lation functions at plane III are detected by a detector 
32 and displayed on a cathode ray tube 33 and/or re- 
corder 34. Once the system has been calibrated, the 
output of the detector 32 provides a direct indication of 
wave height just as in the case of the digital processor of 
FIG. 1(a). There is thus shown an optical processor that 
is directly analogous to the digital processor. 

Although particular embodiments of the invention 
have been described and illustrated herein, it is recog- 
nized that modifications and variations may readily 
occur to those skilled in the art. It is therefore intended 
that the claims be interpreted to cover such modifica- 
tions and variations. 

What is claimed is: 

1. Apparatus for determining heights of ocean waves 
or terrain peaks utilizing synthetic-aperture radar to 
attain a return signal from said waves or peaks, said 
signal having a spectrum of bandwidth f n comprising 

means for compressing said radar signal over defer- 
ent bandwidths, A f n of said spectrum to obtain a 
plurality of intensity images I,(x), I 2 (x), . . . I„(x), 
said bandwidths being centered at frequencies f m , 
// 02 > • • • fan* each spaced by a frequency difference 
of A/o from the first /on and 

means for cross-correlating a selected area of one of 
said images with a corresponding area of at least 
one of the other of said images to determine a 
cross-correlation function value R £ (A/o) for a 
given spacing A/o which is a function of height of 
said waves or peaks. 

2. Apparatus as defined in claim 1 wherein said means 
for cross-correlating said selected area of said one of 
said images with a corresponding area of each of said 
other images includes means for recording a predeter- 
mined component of the resulting cross-correlation 
values as a function of the spacing of the center frequen- 
cies for the bandwidths of the cross-correlated images, 
whereby a measure of the value of a plotted curve for a 
given frequency spacing provides a direct indication of 
height. 

3. Apparatus as defined in claim 2 wherein said com- 
pressing means is comprised of a plurality of filter 
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means of bandwidth A/„ one for each center frequency 
froi> froi* • • • frOn> and each of predetermined width A/„ 
and means for squaring the output signal of each filter- 
ing means to obtain an intensity image signal, and said 
cross-correlating means includes means for cross-cor- 5 
relating a selected area of a first intensity image signal 
IjC*) with a corresponding area of each of said intensity 
image signals to form the cross-correlation values equal 
to fli(x)l k (x)dx for k = 1, 2 thus providing 
cross-correlation values R £ (A/o) for the center fre- 10 
quency spacing, A/o, of said intensity images cross-cor- 
related. 

4. Apparatus as defined in claim 1 wherein said corre- 
sponding means includes optical filtering means operat- 
ing on said radar signal stored on a data film transpar- 15 
ency, and produces a film transparency recording of 
said intensity images, and said cross-correlating means 
is comprised of optical means for superimposing a se- 


lected area of one intensity image on one film transpar- 
ency over corresponding areas of other intensity images 
on another film transparency, and means for detecting 
the cross-correlated values of said images, thus provid- 
ing cross-correlation values R £ (A/o) for the center fre- 
quency spacing, A/o, of said intensity images thus cross- 
correlated. 

5. Apparatus as defined in claim 4 wherein said cross- 
correlating means comprises means for positioning one 
film relative to the other for superimposing one optical 
image at one plane on an optical image in the other film 
in another plane, a source of light for transilluminating 
said films to produce superimposed images, means for 
detecting the product of said images thus superimposed 
at a third plane, and a pupil for restricting light from 
said source to a selected area of said films to produce 
said superimposed images at only said selected area. 
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